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Abstracfc— A larger set of siimiltaneous data froDi the Apollo 12 lunar sur- 
face magnetometer and the Explorer 35 Ames magnetometer are used to con- 
struct a "whole-moon hysteresis curve, from 'vdiich a new value of global lunar 
permeability is determined to be \x> = 1.012 ± 0.006. The correspondxng glo- 
bal induced dipole moment is 2,1 x 10 gauss-cm^ for typical inducing fields 
of lO”^ gauss in the lunar environment. From the permeability measurement, 
lunar free iron abundance is determined to be 2.5 - 2.0 wt. io. Total iron 
ab-undance (sum of iron in the ferromagnetic and paramagnetic states) is 
calculated for two assumed coit 5 )ositional models of the lunar interior: a 

free iron/orthopyroxene lunar con^josition and a free iron/olivine composi- 
tion. The overall lunar total iron abundance is determined to be 9*0 - ^.7 
wt. io^ Other lunar models with a small iron core and with a shallow 
iron-rich layer are discussed in light of the measured global permeability. 
Effects on permeability and iron content calculations due to a possible lunar 
ionosphere are also considered. 


- 1 ' 



IKTBODUCTIOK 


Theoretical calculations of vAiole-moon iron abundance have been made 

by several investigators, often based on meteorite conroositional models. 

Previous estimates for total iron abundance have generally been 10^ by 
> 

weight (Urey, 1962; Reynolds and Summers, 19^9; Urey and MacDonald, 1971; 
Wanke et al., 1973)* 

Enqplacement of a neti/rark of magnetometer experiments on the Imar 
surface by Apollo astronauts in 1969-1972 has allowed investigation of 
lunar iron content and magnetic penaeability, using ^ situ magnetic field 
measurements made simultaneously by instruments on the lunar surface and 
in orbit near the moon. Prom these simultaneous measurements whole-moon 
hysteresis curves can be constructed, from which global magnetic permeability 
can be calculated. (Beharaion (1968) previously had calculated an upper limit 
of 1,8 for global permeability using measurem.ents of a single lunar orbiting 
magnetometer.) The dual-nagnetometer method was first eitroloyed by D;|/’al and 
Parkin (l97l)j who calculated whole -moon relative magnetic permeabi?..:' to 
be = 1.03 i 0.13. The uncertainty in this value was later reduced to 
ji, = 1.029 ± Parkin et al, (1973). This measured permeability was 

used to calculate ferromagnetic free iron in the moon to be 5 - ^ wt.^, 
and total iron in the moon to be 9 - ^ wt.^. 

In this paper the earlier work of Parkin et al. (1973) is extended. 

New values of global permeability, free iron, and total iron are determined 
using more data, improved statistical techniques, and better quantitative 
knowledge of the plasma -magnetic field environment of the moon. Ifew hys- 
teresis curves are plotted using a total of 2703 sets of magnetic field 
averages of data measured simultaneously by the Apollo 12 lunar surface 
nagnetometer and the lunar orbiting Explorer 35 magnetometer. Slopes of 



the hysteresis CTjrves are used to calculate magnetic permeability of the 
moon, from vdiich we calculate free iron abundance in the lunar interior. 

Then we calculate total iron abundance of the moon for assumed free iron/ 
paramagnetic mineral coirpositions of the lunar interior. Also, we discuss 
other lunar models, one with an iron core and another vrith an iron-rich 
layer, in light of the measured global lunar permeability. Finally, Russell 
et al. ( 197 ^) have recently made permeability calculations using data from 
a single magnetometer, the Apollo 15 subsatellite magnetometer orbiting at 
an altitude about 100 km above the moon. The results to date indicate the 
possible existence of a lunar ionosphere betirreen the lunar surface and the 
sub satellite orbit. We consider the effects of such an ionosphere, should 
its existence be proved, upon the lunar magnetic permeability and iron 
abundance results. 


EXEERIlvEKTAl TECKl'KOUE 

During times \ihen the moon is immersed in an external magnetic field 
which is spatially and teuporally taiiform, and plasma interaction effects 
are minimized,' the total magnetic field B meastired at the lunar surface by 
an Apollo magnetometer is expressible as 

B = n H = H + ( 1 ) 

where H is the external magnetizing field and M is the magnetization field 
induced in the permeable lunar material (see Fig. l). The relative magnetic 
permeability is |i = 1 + 4itk, where k is rr.? -'netic susceptibility in emrn/cm^, 
(in equation (l) the remanent field at the surface site is subtracted out 
for siit^jlicity, ) The total surface field B is measured by a lunar surfa'-'r 
magnetometer (LSM) located at the Apollo 12 site on the moon (selenogranhic 
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coordinates 23.4^ W. longitude, 3*0° S, latitude), LSM instrument proper- 
ties are described in detail by Dyal et al. (19T0) and summarized in a 
conpanion paper in this volume (Dyal et al., 197^)* Simultaneous measure- 
ments of the geomagnetic tail field H are made by the lunar orbiting 
Explorer 35 Ames magnetometer, idiich orbits the moon ^ri.th 0.5 periselene 
and 5 Rjj^^ aposelene at a period of 11.5 hours. Characteristics of the 
Explorer 35 mgnetometer are outlined by Sonett et al. (196?) • 

To determine iron content and magnetic permeability of the moon, we 
first construct a lunar B-H hysteresis curve using simultaneous measurements 
pf the magnetizing field H external to the moon and the total magnetic in- 
duction B on the lunar surface. The slope of the hysteresis curve is then 
determined and the viiole-moon initial permeability is calculated. This 

value is called the "initial" permeability because it is associated -vath 

/ -U 

the B-H hysteresis curve at very small magnetizing field H (order of 10 Oe). 
At these small field values, the characteristic "S" shape of the hysteresis 
curve degenerates to a straight line through the origin (Ellwood, 193^). 

Then from the global permeability'’, iron abundance in the lunar interior is 
determined as a function of thermal and compositional models of the lunar 
interior . 


MAGHETIC ENVIE0M.ffiNT IH THE GE0L4AGHETIC T-^iIL 
In different regions of a lunar orbit (see Fig. 2), the magnetic 
en’vlronment of the moon can have distinctly different characteristics. Con- 
ditions desirable for analysis of lunar permeability and iron content exist 
in regions of the geomagnetic -bail where the earth’s magnetic field is 
spatially uniform and teii5>orally constant, and effects of plasma currents 
are negligible. 
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The geomagnetic tail is formed due to an interaction of the earth's 

permanent dipole field with the plasma flowing radially outward from the 

sun at an average velocity of ~ 400 km/sec; in effect, the earth’s field 

is swept back into a cylindrical region (the geomagnetic tail) on the earth’s 
> 

antisolar side. At the distance where the moon’s orbit intersects the tail, 
the field magnitude is ~ 10 gammas (10 Oe), Substructure of the tail 
consists of two "lobes” : the upper or northmrd lobe has its magnetic field 

pointing roughly to^-jard the earth, whereas the lovrer lobe field points avra-y 
from the earth. The moon can pass through either or both lobes, depending 
ttpon the characteristics of the particular orbit, the geomagnetic dipole 
axis orientation, and perturbations of the geomagnetic field by solar Tri-nd 
pressures. 

Data must be carefully chosen in the geotail so that interaction modes 
and induction modes other than global magnetization, are minimized. In 

i 

general, it is possible that the total measured surface field E can include 
field contributions in addition to the geomagnetic field H and the magneti- 
zation field M induced in permeable Ixmar material. These other possible 
contributing field modes, some of which are very isiportant in regions of 
the lunar orbit outside the geomagnetic tail, are the following: (l) the 

steady remanent field due to permanently magnetized subsurface materials 
nc ar the Apollo magnetometer site, (2) the toroidal field ^ resulting from 
transverse magnetic (TM) induction, (-3) the poloidal field B^p resulting from 
transverse electric (TE) induction, (4) the field associated with dynamic 

interaction between plasma flow and the above lunar fields, and (5) the field 
^ associated with plasma diamagnetic currents. We examine here each of 
these possible fields to assess its impact on our study of the lunar magneti- 
zation field. 
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The remanent field ^’d.U not affect the hysteresis slope measurement. 
Since is constant, its components added to equation (l) would, for 
hysteresis curves plotted for the separate field conponents of B and H, 
simply result in a shift of the hysteresis curve aT‘;ay from the origin with- 
out altering the slope. Indeed, the radial component of B at the Apollo 22 

site has been subtracted out of the lunar hysteresis curve presented in 
Fig. 3. 

The toroidal mode B^ would result from currents driven in the lunar 
interior by a motional electric field == V x H, ^■ihere V is relative velocity 
of the moon with respect to the ambient bulk plasma flow. This type of field 
has not been detected within experimental error by lunar magnetometers (Dyal 
and I^rkin, 1971) even in the solar wind, where V ~ kOO km/sec, much faster 
than the moon could move relative to the geomagnetic tail field during quiet 
times (lunar orbital velocity is 1 ka/sec and geomagnetic tall motional 
velocity is generally - 70 km/sec (Milualov et al., 1970)). Furthermore B 
is everyiv-here toroidal to and tangential to the lunar surface (Schubert 
and Schwartz, 1969; Sill and Blank, 1970); therefore, contaminating effects 
of this induction mode can be neglected by use of only the radial components 
of B and H in our hysteresis curve analysis. 

The induced poloidal field B^ results when time-dependent fluctuations 
of the field external to the moon O^St) drive eddy currents in the lunar 
Interior. B^ is the dominant induction mode during times when the moon is 
immersed in the turbulent, rapidly changing solar wind field; however, in 
the limit of low-frequency or small-amplitude driving field fluctuations, 
as in the tail, 3^Bt ** 0 and poloidal induction -vanishes. Therefore data- 
selection restrictions are placed on peak-to-peak variations of surface and 

external fields (to be discussed later) to eliminate data obvlouslv 
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contaminated by eddy current fields. Furthermore, since global induced 
B fields have dipolar symmetry about the direction of BH/dt, rather than 
H, poloidal contamination would tend to scatter hysteresis curve data rather 
than change the hysteresis curve slope. 

When the moon is in the free- streaming solar wind outside the geomag- 
netic tail, the dynamic interaction fields (B^) caused by the interaction 
bet;'reen solar i/ri-nd plasica and the lunar sui'face fields have been found 
to be ^ l6 gamnias at the Apollo 12 site (Dyal et al,, 1973)* Plasma stream- 
ing pressure responsible for compression of surface fields is (l/2)N^mpV^, 
where 17^ is proton number density, mp is the proton mass, and V is the 
plasma b\ilk speed. We assume that the typical plasma speed in the tail is 
characterized by the average magnetopause motions which are about 70 tan/sec 
(Mihalov et al. , 1970). Deep in the magnetotail lobes the plasma density 

is substantially below the density of O.l/cm^ measured in the plasma sheet 

\ 

(Rich et al., 1973 )• Therefore the plasma-field interaction in the geotail 

3 

should be at least 10 times weaker than in the free- streaming solar wind 
and can be neglected. 

Another possible contaminating field is that due to plasma diamagnetism. 
As plasma particles spiral about the magnetic field lines H, their motion 
induces a field (B^^) idiich opposes H. Plasma diamagnetism in the solar 
wind has been measured (Colburn et al. , (19^7)} Ness et al. , (1967)) as a 
field magnitude change of about 1,5 gammas by comparing F^lorer 35 measure- 
ments in the solar wind with those in the plasma void on the antisolar side 
of the moon. Explorer 35 is unable to measure 3^ directly in the plasma 
sheet of the geomagnetic tail, however, since there is no well defined 
plasma void created by the moon in the tail comparable to the antisolar 
cavity formed in the solar wind. .An .Apollo magnetometer on the lunar surface 
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should be shielded (at least partially) from extralunar diamagnetic effects; 
indeed, ve have examined extensive data taken •when the moon is near the 
neutral sheet, and ■we have seen differences as much as 2 gammas in fields 
measured by Apollo 12 LSM and Explorer 35. These field differences are 
directed such that they oppose the direction of the earth’s field H, charac- 
teristic of diamagnetic fields. We ha've found that we can minimize the 
inclusion of plasma diamagnetic fields in our hysteresis curve analysis by 
eliminating all data for ■which the measured Explorer 35 field magnitude 
|h| < 6xl0"^ Oe. This data selection criterion is applied so that only data 
from regions deep in the tail lobes are utilized. 

Diamagnetism in the external environment produces the same effect in 
hysteresis-curve analysis as does paramagnetism in the lunar interior, i.e., 
paramagnetism determined for the moon is only relative to that of the I'unar 
environment. Therefore if diamagnetic field data are included in the radial 
B-H hysteresis curve analysis, the measured global permeability \-7ill be 
higher than the true lunar permeabili'fcy. In our earlier -work (Parkin et al., 
1973 the calculated permeability -^vas higher than in this paper, due in 
part to inclusion of some plasma sheet data in the analysis. 

MAGMETIC PERIfflABILITI OF THE MOON 

Magnetic permeability and iron ab-undance of the noon are calculated 
by analysis of magnetization fields induced in the permeable material of 
■the moon. When the moon is immersed in an external field it is magnetized; 
the induced magnetization is a function of the distribution of permeable 
material in the interior. Under the assumption that the permeable material 
in the moon is predominately free iron and iron-bearing minerals, the lunar 
iron abundance can be calculated from the lunar permeability for assumed 
conpositional models of the interior. Since the 
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amotmt of iro2i present in the Itmar interior should be consistent with the 
measured global magnetic permeability, the permeability in effect places a 
constraint on the physical and chemical composition of the moon’s interior. 

In this section we calculate global magnetic permeability of the moon. Lunar 
iron abundance will be determined from the permeability results in the follow- 


ing section. 

Theory 

For the lunar permeability and iron abxmdance analysis we use a two- 
layer lunar model (see Fig. l). We assume that free iron and chemically 

O , 

combined iron in the Fe (ferrous) state are responsible for the magnetic 
permeability of the lunar interior, and \re model the moon with a hoanogeneous 
paramagnetic rock mtrix (olivine and orthopyroxene models "^-ri-ll be used), 
in which I'ree metallic iron is uniformly distributed. Since the suscepti- 
bility of free iron changes several orders of magnitude at the iron Curie 
teaaiperature (T^,)} a t'vao-layer model is used \rit'a the core-shell boundary 

at the Curie isotherm. For R>R,T<T, and for R ^ R , T > T . There- 

c c c c 

for 

fore^R > R^ any free iron mil be ferromagnetic, while at greater depths 

where T > T , the free iron will be paramagnetic, 
c 

For the two-layer lunar permeability m.odel illustrated in Fig, 1 the 
total field at the outer surface of the sphere is expressed (see Parkin 


C O CL-L. 


*1 rrro ^ # 

1 j / • 


B = H (1 + 2G) X + H (i-G) y + H^(l-G) z 


where 


G - 


(2TVfl)(n5_-l) - x3(n-i)(2u^+l) 

(2r,+l)(u.^+2) , 2 x3(i-1)(ui-1) 


’ (2) 


( 3 ) 



Here t) = and are relative permeability of the shell and core, 

respectively. Again, are coirponents of the geomagnetic tail field 
external to the moon, measured hy the Itmar orbiting Ames Explorer 35 mag- 
netometer; B is measured by the Apollo 2 .^. limar surface magnetometer. The 
> 

permeability exterior to the sphere is |j,^ = 1, that of free space; X = 

R and R are radius of the core and the moon, respectively. Equation (2) 
c m 

expresses the total surface field in a coordinate system which hs-s its origin 
on the lunar surface at an Apollo magnetometer site : x is directed radially 

outvrard from the lunar surface, and y and z are tangential to the surface, 
directed eastward and north^'jard, respectively, 

A plot of any component of equation (2) ttIII result in a B - H hysteresis 
curve; Equation (3) relates the slope of the hysteresis curve to the lunar 
permeability. The average whole-moon permeability |j, is calculated from the 
hysteresis-curve slope by setting = ii^ = P- in equation ( 3 ): 


- iL 


p, + 2 




Hysteresis curve data selection criteria 

The most recent data reduction has resulted in the global lunar 
hysteresis cu3nfe shown in Fig. 3- This curve has been constructed using a 
total of 2703 simiiltaneous Apollo 12-Explorer 35 magnetometer data sets, 
measfirecl. during four orbits of the moon through the geomagnetic tail, 

(in the earlier permeability calculation (Parkin et al. , 1973) a- smaller 
quantity of Z03 data sets was used, with less strict data selection criteria. 

The present hysteresis curve data have been carefully selected to mini- 
mize the aforementioned contaminating induction and interaction field modes. 
For reasons described in the previous section, contamination from plasma 
diamagnetism (B^^) in the exti’alur'ar environment is minimised by eliminating 
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data points for 'which the magnitude of the external magnetizing field 
|Hj < 6x10“^ Oe. 

Since poloidal eddy-cnrrent induction (Bp) is dependent upon time rate 

of change of the external field (BH/Bt), contamination firom the poloidal 
> 

mode is minimized by restricting variations in the driving field H, Two- 

minute inter'vals of H and B data (vri.th |h| > 6x10 ^ Oe) are examined using 

a conputer program designed to select intervals during which the Explorer 

-5 -5 

35 and Apollo 12 data peak-to-peak variations are < 2x10 Oe and < 1x10 Oe, 
respectivelj'-. Data which qualify under these criteria are then averaged 
over the t'wo-ninute intervals, and the averages are plotted to construct 
hysteresis curves. For likely lunar electrical conductivity profiles (Dyal 
et al;, 197 ^) > the peak-to-peak criteria iirposed on da-ta, limit poloidal 
field induction to magnitude <0.7 gamma. We average data over t'vro-minute 
intervals, however, so that the possible error due to poloidal induction 
is much less than 0.7 gamma. Furthermore, since Bp is not a function of 
H but rather of B^Bt, errors due to poloidal induction could add scatter 
to the linear hysteresis curve, but not change its slope. The same data- 
selection process has also been repeated for six-minute intervals of H and 
B data, 'With similar results. These peak-to-peak criteria have been, ecployed 
both to minimize poloidal field contamination and to insure uniformity of 
the field K over a region larger than the distance between the Apollo 12 
and Explorer 35 for all selected data.. 

We note here that similar peak-to-neak criteria could not be applied 
to simultaneous Apollo 15 LSM and Explorer 35 data since by the time Apollo 
15 was deployed years after the Apollo 12 data were taken, the Explorer 
35 magnetometer data were contaminated by a spin tone modulation, of period 

_5 

0.8 min, which caused peak-to-peak oscillations as imich as 5x10 ■ Oe 
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during quiet times in the geomagnetic tail. Therefore as yet we have been 

unable to update the Apollo 15 results which -vtc reported earlier (Parkin 

et al., 1973 )j and ve consider only Apollo 12 data in this paper. 

In our calculations of the whole moon permeability we use only results 

from the 'radial components of B and H data for several reasons. When the 

moon is in the geotail the Apollo 12 LSM location on the lunar surface is 

such that on the average, the geotail field has a radial congsonent at least 

twice as large as either y or z tangential component, i.e., H 2 H . 

X y , * 

Also, for the dipolar magnetization field M, M /H = 2 M „/H „ (see 

equations (l) and (2)). Therefore M ^ 4 M , i.e., the resolution of the 

X y, z 

magnetization field is at least four times greater using the radial com- 
ponent. Furthermore two contaminating field modes are minimal in the radial 
component. Toroidal fields are not present in the x-component as explained 
in a previous section, and plasma interaction effects (measured in the solar 
•vrind) are much lov/’er in the x direction at the Apollo 12 site; /B_, 0.2 

and ~ O.T (l»yal et al., 1972). 

Global magnetic perrr.eability and induced dipole moment 

The global lunar hysteresis curve wiiich meets data selection criteria 
discussed in the previous section is shovm in Fig. 3* Apollo 12 averages 
of radial (x) components are plotted on the abscissa with simultaneous 
Explorer 35 averages on the ordinate; 2703 sets of selected two-minute 
averages are used. The gap in the center of the curve is a result of the 
data selection criterion which eliminates data measured when the moon is 
near the neutral sheet. The familiar "S" shape of the hysteresis curve 
degenerates, at these low field values to a straight line (Ellwood, 1934) 
intersecting the origin. The x-component of the remanent field at the 
Apollo 12 site has been subtracted from the Apollo 12 data.. To this data . 
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set a least squares line has been fitted. The best estiraate of the slope 
is 1.008 ± 0.004. 

The least-squares result is obtained as follows. Since both Exq)lorer 
35 and Apollo 12 data have errors assumed to be mutually independent and 


Gaussian in nature, neither data set can be considered to be the independent 
set. For this reason ti-ro linear regression analyses have been used, one 
with B as the independent variable, and the other using H as the indepen- 
dent variable. The ttro regression coefficients have been calculated to be 
1.005 - 0,002 and 1.011 ± 0.002, respectively; the ± 0.002 error limits are 
standard errors. To get the final value we have bisected the angle betvreen 
the two least-squares lines and calculated the slope of the bisector to be 
1.008< The 95^ confidence interval of the regression analyses is ± 0.004; 
therefore our calculated best estimate of the slope is 1.008 ± 0.004, 

Scatter in the Apollo 12 -Explorer 35 data points of Fig. 3 is primarily 
a result of magnetic inhonogeneities between the moon and Explorer 35, small 
contributions .from eddy current fields, and instrumental noise and offset 
errors in the Apollo and Explorer magnetometers. These error soirrces may 
introduce small random fluctuations into the data vrhich will not substan- 
tially affect the slope or intercept of the least-squares line. 

From the slope we calculate global relative magnetic permeability of 
the moon to be y, = 1,012 ± 0,006 using equations (2) and (4). Both extrema 
are greater than 1,0, ir^ilying that the moon, as a vdiole, acts as a para- 
magnetic Of weakly ferromagnetic sphere. 

The global induced dipole moment can be calculated from the global 
permeability. The induced dipole moment is expressed as GE^ H, in units of 
gauss-cm'^; the hysteresis curve slope = 2G + 1. Using = 1.74x10 cm and 


georngnetic tail field H = 10 ' Qe, 


the induced dipole moment is determined 


- 13 - 



IS 3 

to be 2.1x10 gauss-cm . 


IRON ABUnOMCE IK KE MOON 

In this section we calcvilate firee iron and total iron abundances in 
the moon >frora the global permeability result of the preceding section. 
Calculations are in general dependent upon thermal and coirpositional models 
of the lunar interior, and are constrained by the known lunar density and 
moment of inertia. 

Theory 

We refer again to the two- layer model of the moon shown in Fig, 1, 

The core-crust boundary at R = is the iron Curie isotherm (assumed to 
be spherically symmetric); its location is a function of the thermal pro- 
file of the lunar interior. In our calculations xre assume the moon is 
coDposed of a homogeneous mineral (orthopyroxene or olivine) of lunar den- 
sity 3 *3^ gm/cm^. Free iron grains are dispersed uniformly throughout the 
lunar sphere. For R > R^ any free iron is ferromagnetic while at greater 
depths vrtiere T > T^, the free iron is paramagnetic. 

In the outer shell there are both ferromagnetic and paramagnetic con- 
tributions to the total magnetic permeability = 1 + The suscepti- 
bility of the shell is where is "apparent” ferromagnetic 

susceptibility and kj^^ is paramagnetic susceptibility. The ferromagnetic 
con 5 >onent is metallic free iron, assumed to be composed of multidomain, 
noninteracting grains; the paramagnetic component is Fe^’’’ combined in the 
orthopyroxene or olivine rock matrix. The expected pressures and tempera- 
tures in the outer shell are such that the ferroioagnetic susceptibility of 
iron will not be subs'fcantially altered (Eozorth, 1951; Kapitsa, 1955). 

The measured ferromagnetic sxisceptibility of the shell material is an 



apparent value \^ich differs from the intrinsic ferromagnetic susceptibility 
of the iron because of self-demagnetization of the iron grains and the volume 
fraction of iron in the shell. The apparent ferromagnetic susceptibility 
k is related to the intrinsic susceptibility according to 


O.Olp 

^a "" 


(5) 


where N is the shape demagnetization factor of the grains and p is the volume 
percentage of free iron in the lunar material. For spherical iron grains 
N = 4it/3 theoretically, but experimentally this value is found to range from 
3 to 4 (Nagata, I961). We shall use M = 3.5 in our calculations. 

Unlike the free iron, the paramagnetic iron in the rock matrix has a 
continuous susceptibility across the Curie isotherm. For the combined iron 
the susceptibility varies according to the Langevin relation; 

k^^ = nm^/3KT (6) 


where K is the Boltzmann constant, T is absolute temperature, n is the number 
of ions per gram, and m is the atomic moment, m is of the order of a few 
Bohr magnetons e.g., for the Fe ion ra = 5.25 to 5.53 (Nagata, 


1961). 

i'or K ^ the lunar maueriLal is paramagnetj-C only, wxth susceptibxlxty 
= k^^ + ^2c contribution of paramagnetic chemically combined 

iron and is the contribution of free paramagnetic iron above the Curie 
temperature. Again, kg^^ is only an apparent value and is related to kg^, 
the paramgnetic susceptibility 'of free iron, by an equation similar to (5); 
and is dependent on temperature with a relationship analogous to equation 


( 6 ). 
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Compositional and thermal models 

Lunar iron abundance is determined for two con^ositional models of the 

lunar interior. Kecently Kaula et al. ( 197 ^) have determined the lunar 

moment of inertia I to be = 0.3952 * 0.0045, where M and are 

lunar mass and radius, respectively. This value is approximately that of 

2 

a homogeneous sphere of constant density = 0.400); therefore we assume 

a moon of uniform density p = 3.34 g/cm , i.e,, we will use homogeneous 
compositional models. 

In the first conroositional model we consider the Ixmar sphere to be 

composed of orthopjnroxene (xPeSiO^ * (l-x)l/lgSi 02 ) , where x is the mole 

2 + 

fraction of the Fe phase present. Free iron grains are dispersed uniformly 
throughout the sphere. In a second model we replace the orthopyroxene with 
olivine (xFe^SiOj^ • (l-x)l!gg3i0]i^). Pyroxenes and olivines have been reported 
to be major mineral components of the lunar surface fines and I'ock samples 
(Nagata et al,, 1971; Zussman, 1972; Weeks, 1972), with combined iron present 
as the paramagnetic Fe^”*” ion. The ferromagnetic component of lunar sairsples 
is primrily metallic iron which is sometimes alloyed with small amounts of 
nickel and cobalt (Nagata et al. , 1972; Pierce et al., 1971 )• This free 
iron is thought to be native to the moon (because of its low nickel content) 
rather than meteoritic in origin (Strangivay et al. , 1973). Orthopyroxene 
and olivine models are consistent •^•o.th geochemical studies (Urey et al. , 

197 I; Wood et al., 1970; Ringwood and Essene, 1970; Green et al., 1971) aJifi 
geophysical studies (Toksoa, 1974). 

Since the susceptibility of free iron changes several orders of magni- 
tude at the iron Curie tempera tirre (t^), a two-layer model has been used, 
with the core-shell boundary R^ at the Curie isotherm (see FiguD v 10). 
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The Curie isotherm location is determined from the thermal profile used 
for a particular model. Fig. 4 shows the dependence of on depth in the 

moon for hydrostatic equilihrium, superimposed on temperature profiles 
proposed for the moon by several authors; the increased pressure of the 
interior will decrease the iron Curie point by about 7^10 C per atmosphere 
increase in pressure (Bozorth, 1951). For the calculations that follow, vre 
have constructed three temperature models to span the range of tecpierature 
profiles in Fig, 4. Presently we use tzro-layer teirperat'ore profiles for 
simplicity. For model profile the Curie isotherm is spherically sym- 
metric and located at R /R =0,9. Shell and core temneratures are 600 °C 
and l400 °C, respectively. For the model profile the shell is 500 ^C, 
and the core is 1300 while the Curie isotherm boundary is at = 0,85. 

Tenperatures are 300 °C and 70O °C for shell and core of model profile T^, 

which has R /R = 0,7. 
d m 

Normally the ferromagnetism of free iron is -dependent on pressure as 

well as tenperature. At very low magnetizing field strengths such as those 

of the geomagnetic tail, hovrever, the susceptibility of iron is not strongly 

dependent on temperature below the Cxarie point (Bozorth, 1951). Uniajcial 

stress on iron changes its susceptibility (Kem, I961); however, hydrostatic 

stress should not affect the susceptibility (Kapitsa, 1955) unless, at very 

high pressures, there is a change in volume (Breiner, I967). Therefore we 

assume that the susceptibility of uncombined lunar iron is independent of 

pressure. In the outer shell where T < T , we define intrinsic ferromagnetic 

c 

o 

susceptibility of the free iron to be = 12 emu/cm^ (Bozorth, 1951 ), and 
in the core where T > T^, free iron intrinsic paramagnetic susceptibility 
is = 2 X 10~ emu/cm”’ (Berko^rttz and --.leller, 19^9; Bozorth, 1951). 



I 

Global limar iron abtmdance 

Using the information described in previous sections we have generated 
the curves shovm in Pig. 5j wiiich relate free iron abundance (q) and total 
iron abundance (Q) to hysteresis-ciorve slope, as follows. The apparent 
ferromagnetic susceptibility of free iron in the shell and the apparent 

paramagnetic susceptibility of free iron in the core each calcu- 

lated as a function of the free iron abundance (q) using equation (5). The 

shell-core boundary is defined by the Curie isotherm used for a particular 

2+ 

temperature model. In addition, the mole fraction of the Fe phase in para- 
magnetic rock used in the model (orthopyroxene or olivine) is constrained 
by q and the bulk density. Mole fraction of Fe is related to suscepti- 
bility using the experimental data of Nagata et al. (1957 ) olivine and 

Akimoto et al. (1958) for orthopyroxene. Furthermore, using the Langevin 

2-l 

relation (equation (6)), v;e relate susceptibility of Fe ' to temperature 
and thus find Fe^"*" susceptibility in the shell "the core (k^^) foj^ 

our three tenperature models. Then we combine susceptibilities to obtain 
total shell and core susceptibilities, k^ = k^ + k^^^ and ~ + k2c’ 

as a function of q, after which we use k^ and k^ in equation (3) to relate 
free iron abundance q directly to G as sho^m in Fig, 5» G is related to 
the slope of the hysteresis curvet slope = 2G + i. 

We complete Fig. 5 by determining Q, the total iron abundance in the 
moon, as a function of G. Total iron Q is the combined abundance of free 
iron and ferrous iron. From the previously determined mole fraction of 
the phase in the paramagnetic mineral, vre calculate the mass of ferrous 

iron as a function of q. Then we add the ferrous iron mass to free iron 
mass and get the total mass of iron in the moon, and thereby the total iron 
abunda.nce Q, 
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We note that for q > 1 wt. ^ the susceptibility of the moon is dominated 

by the ferromagnetic iron (generally and. therefore 

the relationship between G and q is independent of the coirposition of the 

ferrous coriponent (olivine or orthopyroxene). Total iron abundances differ 
> 

for the tvro compositional models because of the different iron contents of 
olivine and orthopyroxene . 

From our measured hysteresis slope of 1.008 ± O.OOlf we calculate 
G = 0.004 ± 0.002. Using this range of G, we find from Fig. 5, free iron 


aboundance q and total iron abundance Q for each thermal and cocpositional 
model. These iron abundances are summarized in table 1. 

Fig, 5 shows the ranges of free iron and total iron abundances, which 



are functions of teaperatiire in the lunar Interior and are bounded by errors 
in the hysteresis curve slope. The lunar free iron abundance ranges between 
4.5 and 0.5 wt. these limits correspond to thermal profiles Tj^ and T^, 
respectively. These therifial proiiles have been selected to be unper and 
lovrer limiting cases, as can be seen upon comparison of Fig. 4 with thermal 
profile descriptions in the text. We accordingly calculate our free iron 
abundance best value as q = 2.5 ± 2.0 v/t. io. 


Total iron abundance (Q) in the moon is, in addition, dependent upon 
compositional model. For the free iron/orthopyroxene model, upper and lov;er 
limits on Q are 13.? and 11.8 -tvt. respectively; for free iron/olivine, 
limits are 6.8 and 4,3 wt. %, Assumng that the moon is composed of one or 
a combination of these minerals and has a total iron abundance betireen 13 .? 
and 4,3 "wt. we calculate the total iron abundance best value as 
Q = 9-0 ± 4.7 wt. Free iron and total iron abundances are sho^m in 



6 
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We note here that the susceptihilities of both olivine and orthopyroxene 

are about an order of magnitude too snail to account for the measured lunar 

permeability without some ferromagnetic iron present. Also, we can calcix- 

late the tp-T nimum free iron abundance in the moon consistent with the 

hysteresis -curve measurements. To do this, we consider the extreme case 

where the measured ■tJhole-moon permeability is assumed to correspond entirely 

to ferromagnetic iron in the outer shell of the moon where the tenperature 

is below the Curie point. For this case the bulk lunar iron abundance 

0. is 0.9 1 0.5 wt, ^ 0 . This result is independent of Curie-point depth of 

extreme 

our three models, The/lovrer limit placed on the lijnar free iron abundance 
by our analysis, therefore, is O.i)- wt. °!o. 

Considerations of an Iron core and iron-rich layer 

The ■sdiole-moon permeability has also been used to investigate the 
magnetic effects of a hypothetical iron core in the moon. Density and 
moment of inertia measurements for the moon limit the size of such a core 
to less than 500 km in radius (Toks’dz, 1974). If this hypothetical iron 
core were entirely paramagnetic and the surrounding core vrere orthopyroxene 
of average tecperature 1100 "^C the global permeability would be I.OOO 3 . , 
This value is small coirpared to the measured permeability of 1.012 ± 0 . 006 , 
irplying that if such a small paramagnetic iron core exists, its magnetiza- 
tion is masked by magnetic material lying nearer to the siirface. Therefore 
the hysteresis measurements can neither confirm nor rule out the existence 
of a smalJ iron core in the moon. 

An iron-rich layer in the moon has been considered by several investi- 
gators (e.g., Wood et al, , 1970; Urey et al. , 1971; Oast and Giuli, 1972). 

It is possible that early melting and' subsequent differention of the cuter 
several hundred kilometers of the moon may have resulted in the foriration 
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of a high-density, iron-rich layer beneath a low-density, iron-depleted 
crust. Constraints have been placed on an iron-rich layer by Cast and 
Gi\ili (1972) using geochemical and geophysical data (for example, measure- 
ments of lunar moments of inertia). One set of their models consists of 
% 

high-density layers between depths of 100 km and 300 km. At a depth of 
100 km the allowed layer thickness is 12 km; the thickness increases vrith 
increasing depth, to 50 km at 300 km depth. Also presented are a set of 
layers at 500 1cm depth. By using exactly the same considerations as were 
used in the iron abundance calculations, we calculate whole-moon permeabili- 
ties which would be expected from limar models cvlth these iron-rich layers. 
The calculations indicate that iron-rich layers allo>red by geophysical 
constraints as outlined by Cast and Giuli, if wholly above the iron Curie 
temperature and therefore paramagnetic, would yield global permeabilities 
of about 1,00006. As for the case of a small lunar iron core, the magneti- 
zation field of such paramagnetic layers would be masked by ferromagnetic 
materials elsewhere in the moon, and the hysteresis ciorve measurements can 
neither confirm nor rxile out these layers. This conclusion vrould particu- 
larly apply to the Gast-Giuli layers at 500 km depths, \rdiich vrould be almost 
certainly paramagnetic. 

If the iron-rich layers existed shallow enough in the moon to be below 
the Curie tempera tiore and were therefore ferromagnetic, then the measured 
global permeability would be about 3.5. This value is well above the upper 
limit for the actmal measured permeability of 1.012 ± O.OO6, and therefore 
the Gast-Giuli layers can be ruled out if they are cool enough to be ferro- 
magnetic. It is important to note that the high-density layers discussed 
by Cast and Giuli (l9?2) can be thought of as limiting cases and that there 
are innuR'.erable less dense ana thinner layers ^-hieh are allovred by geophysica 
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geochemical and magnetic constraints. 

Effects of a possible l\mar ionosphere on permeability and iron abundance 
calculations 

Measjirements of the Apollo 15 sub satellite magnetometer in the geomag- 
netic tail have indicated the possibility that an ionosphere exists in the 
region betvreen the lunar surface and the subsatellite 100 Ion mean altitude 
(Russell et al., 197^). Also, results from the Rice University supra thermal 
ion detector experiments indicate that chazged particles measured on the lunar 
surface could be from a lunar ionosphere (Lindeman et al., 1973 )• If ^ 
global lunar ionosphere does exist, then it could form a diamagnetic region 

around the moon which could lovzer the geotail field K to a smaller magne- 

♦ 

tizing field H’ and in turn result in a smaller total stirface field B’ . 

Then plots of Apollo 22 data versus Explorer 35 “lata would be plots of 
B • vs H , causing the measured permeability to be lower than the true 
lunar permeability (in contrast to the plasma diamagnetism effect discussed 
earlier, z-/hich would tend to make the measured permeability higher). 

If the assumption of Russell et al. (l97^)j that a spherically 
symmetric diaiDagnetic ionosphere fills the entire region between the lunar 
surface and the subsatellite 100 km altitude, were T?ound to be correct, 
then our global lunar permeability irould be adjusted up^-ra.rd. We have 
determined the adjusted permeability and iron abuiidance values by modifying 
equation (,3 j for a two-layer (shell/core) permeability model, where in this 
case the shell represents a lunar ionosphere of 100 km thickness and homo- 
geneous permeability and the core represents a lunar sphere of bulk 
permeability |i. The modified version of equation (3) is 
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( 7 ) 


g = = 

(2n^+fj,)(fJ,^+2) - 2X^(fi^-|j,)(pL^-l) 

■where X = K /E ; R and R are radius of the core (here the lunar glohe) 

-> c' s c s 

and the outer radius of the shell (ionosphere), respectively. From the 
results of-Russell et al, (197^) we obtain g = - 0.0101 ± O.OO 39 . Then 
using soltitions of equations ( 3 ) and ( 7 ) find that our permeability best 

■value ■would be adjusted up'^’ra.rd slightly from 1.012 to I.OI 7 , provided a 
lunar ionosphere exists. The corresponding free iron abundance best ■value 
■would be adjusted upv?ard from 2.5 to 3*9 Total iron content best 

values would be adjusted do\ra^S'ra,rd slightly from 12.8 to 12.1 ■vrt. ^ for the 
free iron/orthopyroxene model, and from 5.5 to 4.7 wt. ia for the free iron/ 
oli^vine model. Be’bailed calculations of the effects of a lunar ionosphere 
on lunar iron content determinations vrill be deferred until a more detailed 
ionospheric model is presented. 
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SUMMARY AND CONCLUSIONS 


(1) Simultaneous measurements by lunar magnetometers on the surface 
of the moon and in orbit around the moon are used to construct a whole'-moon 
hysteresis ^'urve, from which the global lunar relative magnetic permeability is 
determined to be 1 .012 ± 0 *006 . 

(2) The corresponding global induced magnetization dipole moment is 

77 

expressed a = 2.1 x 10 H. For typical geomagnetic tail fields of H = 10 Oe, 

18 3 

the corresponding induced dipole moment is 2.1 x 10 gauss-cm . 

(3) Both error limits on the magnetic permeability value are greater than 
1.0, implying that the moon as a whole is paramagnetic and/or weakly ferromag- 
netic . Assuming that the ferromagnetic component is free metallic iron of multi- 
domain, noninteracting grains, the free iron abundance in the moon is calculated 
to be 2.5 ± 2.0 wt. %. 

(4) A free iron abundance extreme lower limit of 0.4 wt. % is calculated 
under the assumption that the global susceptibility is due entirely to free iron in 
the ferromagnetic state. This lower limit is independent of composition of the rock 
matrix making up the bulk of the moon. 

(5) Total iron abundance in the moon is determined by combining free 
iron and paramagnetic iron components for two assumed lunar compositional models , 
of orthopyroxene and olivine. For an orthopyroxene moon of overall density 

3,34 g/cm^ with free iron dispersed uniformly throughout the lunar interior, 
the total iron abundance is 12.8 ± 1.0 wt. %. For a free iron/olivine moon the 
total iron abundance is 5.5 ± 1.2 wt. %. A summary of iron abundance calcula- 
^tions is given in Table 1. Using extreme upper and lower limits in Table 1, the 
overall total iron abundance is expressible as 9.0 ±4.7 wt. %. 

(6) Lunar models with a small iron core and with an iron-rich layer are 
also discussed using the measured global lunar permeability as a constraint. A 
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small pure iron core of 500 km radius (the maximum size allowed by lunar 

density and moment of inertia measurements) , which is hotter than the iron Curie 

\ 

point (T>T^) , would not be resolvable from the data since its magnetization field 
would be small compared to the induced field we measure. Similarly, an iron- 
rich layer in the moon could not be resolved if the iron is paramagnetic, i.e. , the 
iron is above the iron Curie temperature. Cast and Giuli (1972) have proposed 
a family of high-density layer models for the moon which are geochemically 
feasible. If these models are ii-on-rich layers lying near the lunar surface so 
that T<T^, the ferromagnetic layers would yield a global permeability value well 
above our measured upper limit. Therefore we conclvide that such shallow 
iron-rich-layer models are not consistent with our magnetic permeability 
measurements . 
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TABLE 1 


IRON ABUNDANCE OF TIIE MOON AS A FUNCTION 
OF THERMAL AND COMPOSITIONAL MODELS 



Free Iron Abundance, wt. 

Total Iron Abundance, wt. % 

•/ 

thermal 

mode]. 

c ompo s i t 

model 

L 

^2 


^3 

^2 


orthopyroxene 

1.0 ± 0.5 

2.0 ± 1.0 

3.0 ± 1.5 

13.4 ± 0.3 

13.0 ± 0.5 

12.6 ± 0.6 

oI ivine 

1.0 ± 0.5 

1 

2.0 ± 1.0 

3.0 ± 1.5 

6.5 ± 0.3 

5.9 ± 0.7 

5.3 ± 1.0 


FIGURE CAPTIONS 


Fig. 1 Magnetization induction in the moon. When the moon is immersed in a 

uniform external field H (in this case the steady geomagnetic tail field) , 
a dipolar magnetization field M is induced in permeable material in the 
lunar interior, with the dipole axis of M aligned along the direction of 
H. The total magnetic field near the moon is B = H + 4 7t M . The magnetic 
permeabilities of the two layers are and p 2 » and for regions outside 
the moon, p = = 1 (free space) . H is measured by the lunar orbiting 

Explorer 35, whereas B is measured by an Apollo lunar surface magneto- 
meter (LSM) . Measurements of B and H allow construction of a B-H 
hysteresis curve for the sphere, from which permeability and iron 
. abundance can be calculated. 

Fig. 2 Magnetic environment of the moon during a lunar orbit, with emphasis 

on the geomagnetic tail region. The plane of the lunar orbit very nearly 
coincides with the ecliptic plane of the earth's orbit. The earth's 
permanent dipole field is swept back into a cylindrical region known as 
the geomagnetic tail; at the lunar distance the field magnitude is ~10 
gammas (10 ^ oersteds). Also shown is the shock surface ("bow shock") 
created by supersonic flow of solar wind plasma past the geomagnetic 
tail. Substructure of the tail consists of tv/o "lobes"; the upper or 
northward lobe has its magnetic field (dark arrow) pointing roughly 
toward the earth, W'hereas the lower lobe field points away from the earth. 
I'he moon is immersed in the tail about four days of each orbit; the moon 
can pass through either or both lobes, depending upon the characteris- 
tics of the particular orbit and the orientation of the earth's magnetic 
dipole axis . 
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Fig. 3 


Fig. 4 


Fig. 5 


Hysteresis curve for the moon. Data points are 2703 simiiltaneous 
2~minute averages of radial components of the external geomagnetic 
field data H (measured by the lunar orbiting Explorer 35 Ames magneto- 
meter) and total magnetic induction pH (B is measured by the Apollo 
12 lunar surface magnetometer) . Data points are selected from four 
lunations of measurements made when the moon is immersed in the uniform 

geomagnetic tail field, far from the neutral sheet in the tail. In this 

-4 

low-external-field regime (~ 10 gammas or 10 Oe), the hysteresis curve 
is linear and is fitted by a least-squares line of slope 1.008 ± 0.004. This 
slope corresponds to a whole-moon magnetic permeability of 1.012 ± 0.006. 
The least-squares line intersects the origin exactly in this figure because 
, the vertical-axis intercept (the radial or x-component of the remanent 
field at the Apollo 12 site) was subti'acted out after the least-squares fit 
was made . 

Temperature profiles for the lunar interior published by various authors: 
(1) Hanks and Anderson (1972); (2) Dyal et al. (1974); (3) Toksdz 
(1974); (4) Sonett et al. (1971). Supei-imposed is the pressure- 
dependent Curie temperature for iron versus depth in the moon. 

Iron abundance in the moon as a function of global hysteresis curve 
slope. Free iron abundance (q) and total iron abundance (Q) versus 
the parameter G (hysteresis curve slope equals 2G+1) for three 
temperature profiles described in the text. Total iron abundance is 
shown for two lunar composition models, orthopyroxene and olivine. 
Arrows below the horizontal axis show the range of the parameter G, 
experimentally determined from' the hysteresis curve slope: G = 0.004 
± 0.002. The shaded region defines the allowed values of free and total 
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iron abundances , bounded by hysteresis curve error limits and by 
thermal models T| and T^. 

Fig . 6 Summary of global lunar magnetic permeability and iron abundances . 
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LUNAR IRON ABUNDANCE VS HYSTERESIS SLOPE 
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